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We address the issue of the reducibility of the dynamics on a multilayer network to an equivalent 
process on an aggregated single-layer network. As a typical example of models for opinion formation 
in social networks, we implement the voter model on a two-layer multiplex network, and we study 
its dynamics as a function of two control parameters, namely the fraction of edges simultaneously 
existing in both layers of the network (edge overlap), and the fraction of nodes participating in 
both layers (interlayer connectivity or degree of multiplexity). We compute the asymptotic value 
of the number of active links (interface density) in the thermodynamic limit, and the time to reach 
an absorbing state for finite systems, and we compare the numerical results with the analytical 
predictions on equivalent single-layer networks obtained through various possible aggregation pro¬ 
cedures. We find a large region of parameters where the interface density of large multiplexes gives 
systematic deviations from that of the aggregates. We show that neither of the standard aggregation 
procedures is able to capture the highly nonlinear increase in the lifetime of a finite size multiplex at 
small interlayer connectivity. These results indicate that multiplexity should be appropriately taken 
into account when studying voter model dynamics, and that, in general, single-layer approximations 
might be not accurate enough to properly understand processes occurring on multiplex networks, 
since they might flatten out relevant dynamical details. 

PACS numbers: 


I. INTRODUCTION 

Real-world interactions often happen at different levels 
and are therefore properly modelled by means of multi¬ 
layer networks. Such multilayer approaches m have 
been applied to fields ranging from energy infrastructure 
0 and transport 00, to epidemiology [8]. The multi¬ 
layer set up can either describe interconnected networks 
with nodes of the same nature in each layer, but inter¬ 
acting with nodes of different nature in a different layer, 
or a multiplex structure with nodes of the same nature 
interacting via a different network in each layer. In any 
case, a central methodological question is that of multi¬ 
layer reducibility, that is, when the multilayer framework 
is really needed to explain new phenomena, or when the 
system description can be reduced to an appropriately 
aggregated or reformulated single-layer network. An in¬ 
teresting contribution in this direction has recently come 
from the study of the structural reducibility of multilayer 
networks, i.e. of the possibility of aggregating some of 
the layers of a multi-dimensional network while preserv¬ 
ing its distinguishability from the corresponding single¬ 
layer aggregated graph [9]. Although some recent works 
have pointed out that multiplex dynamics can be intrin¬ 
sically different from their equivalent single-layer coun¬ 
terparts puns, little attention has been devoted so far 
to the problem of reducing a process taking place on a 
multilayer network to a dynamically equivalent process 
on an appropriate single-layer network. In this paper we 
address this general question, by implementing the voter 
model on a multiplex networks and by studying the re¬ 
ducibility of its dynamics as a test case. 

The voter model m is a nonequilibrium lattice 
model [TBj which gives a standard framework for studying 


the influence of social imitation on the process of opinion 
formation im. a basic question considered in this con¬ 
text is when and how the system reaches an absorbing 
state with all the interacting nodes in the same state, or 
when an active dynamical situation of coexistence of dif¬ 
ferent states prevails. The answer to these questions is 
known to depend crucially on the network structure and 
on the update rules employed p°8d22] . The voter model 
has also been instrumental to understand fundamental 
phenomena in coevolution dynamics in which node states 
and network structure have coupled dynamical evolution 
with two different time scales [2311271 . In terms of com¬ 
parison with real data, a metapopulation voter model 
has been recently shown to be able to account for voting 
patterns in the US general elections [28] . 

Quite frequently, social interactions happen on differ¬ 
ent concurrent contexts, so that any model of social im¬ 
itation based on a single-layer representation of social 
relationships should probably be regarded only as a first- 
order approximation of a potentially more complex dy¬ 
namics. In order to better mimic the multidimensional 
nature of social interactions, we consider here a multiplex 
voter model, where agents interact at two distinct layers. 
We assume also that the system exhibits a certain level 
of multiplexity, meaning that a fraction q of the agents 
is present on both layers [26] . As a consequence, any 
change in the state of those agents on either of the two 
layers is automatically transferred to the other layer, ef¬ 
fectively coupling the voter model dynamics taking place 
on each of the two networks. By taking into account 
heterogeneity in the participation of agents to layers we 
aim at reproducing an interesting feature of real-world 
multilevel social systems [29] . In particular, it has been 
shown that the percentage of users of an online social 
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FIG. 1: (Color online) In the multiplex voter model each node i on layer a is associated to a state s(°^(t) = =L1 (the values +1 
and —1 are respectively indicated in the figure by red and blue), which evolves according to one simple rule: select one of your 
neighbours x on layer a uniformly at random, and copy its state, i.e. set s\t*\t + 1) = s\?\t). For instance, the state of node 
l on layer 2 changes from = — 1 (blue, panel a) to s =s +1 (red, panel b) since l selected its neighbour m on layer 2 and 
copied its state on that layer. The presence of inter-layer edges (dashed lines) denotes identification between nodes at the two 
layers. An example is the inter-layer edge connecting the two replicas of node i. In this case, any change in the state of node 
i at one layer will enforce a change of its state on the other layer. In the figure, the state of node i flips from red (panel a) to 
blue (panel b) on both layers, as a consequence of the interaction of node i with its neighbour j on layer 1. 


network (e.g., Twitter) which has also an account on an¬ 
other online social network (e.g., Facebook or Linkedln) 
lies somewhere between 30% and 70%, as reported by the 
social media matrix periodically published by the Pew 
Research Center [50]. These empirical findings confirm 
the necessity to consider systems whose multiplexity is 
neither 0 nor 1, but somewhere in between. A second im¬ 
portant feature of real-world social interactions is the fact 
that networks corresponding to different contexts share a 
number of common links, as can be found for instance in 
data from online games: analysis of the PARDUS online 
society in mi EH yields significant overlap between des¬ 
ignation of other players as friends, correspondents, or 
trading partners. This property of multilayer networks 
can be quantified by the edge-overlap parameter [33]. 

In this work we address a fundamental question, that is 
whether the coupling of several voter models into a multi¬ 
plex dynamics, motivated by the multiplicity of contexts 
influencing real interactions among individuals, gives rise 
to qualitatively different phenomena or effects than those 
observed in the classical single-layer setup. In particular, 
we want to establish whether multiplexity makes any dif¬ 
ference for a simple social dynamics like the voter model, 
or if instead the multiplex voter model can indeed be 
reduced to an equivalent voter model dynamics on an 
appropriately constructed single-layer network. 

A preliminary result on the question of the multi¬ 
layer reducibility exists for a bilayer of uncorrelated net¬ 
works evolving according to the voter model with adap¬ 
tive links [26] . but in the limit of zero edge overlap. In 
this model each layer is associated with a network plas¬ 


ticity parameter, that controls the rate at which relations 
among agents are rearranged, such that if the values of 
plasticity at the two layers are sufficiently different the 
system displays a network shattered fragmentation. On 
the other hand when the layers have the same plasticity 
one finds that the results coming from a pair approxima¬ 
tion in the thermodynamic limit are equivalent to those 
of an appropriate aggregated single-layer network [26]. 
We build upon this, in time independent networks, ex¬ 
amining both the asymptotic properties of the thermo¬ 
dynamic limit and the characteristic times to reach an 
absorbing state for finite systems. We also consider dif¬ 
ferent schemes to obtain a possible equivalent aggregated 
single-layer network and highlight the importance of the 
degree of multiplexity and edge overlapping. 

The paper is organized as follows. Sect II introduces 
our voter dynamics multilayer model. Sect. Ill describes 
our numerical findings. These are compared in Sect. IV 
with single-layer theoretical results for two natural meth¬ 
ods of collapsing the multilayer into an aggregated single¬ 
layer network. Sect. V analyzes an optimum single-layer 
reduction method. Conclusions are discussed in Sect. VI. 


II. THE MODEL 

We consider here the case of a 2-layer undirected un¬ 
weighted multiplex network, described by the pair of bi¬ 
nary adjacency matrices {A^\ where = {a^}, 

and = 1 if and only if node i and node j are con¬ 
nected by a link at layer cp and zero otherwise. On each 
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layer we have N nodes. A parameter of interest in this 
study is the average edge overlap uj, that is the probabil¬ 
ity that an edge is present on both layers: 


LjJ 


y a [1] a [2] 

2^jj a i 3 a i 0 

2 K 


(i) 


where 


k = EE 




j>i 1 + a ij a ?j 


is the number of edges of the graph obtained by aggre¬ 
gating the two layers into a single one [33 . Notice that 
uj = 0 only if each edge exists in exactly one of the two 
layers, but not in the other one, while we have u = 1 
only if all the edges exist on both layers. 

Each node i on layer a is associated to a binary state 
variable where s^\t) can be either +1 or — 1. 

Moreover, we assume that a fraction q of the N nodes 
participates in both layers, requiring that if i is one of 
these qN nodes then its state at the two layers will be 
identical at the end of every update. Such nodes partic¬ 
ipating in both layers are chosen randomly at initializa¬ 
tion. We can think of the parameter q as the interlayer 
connectivity or the degree of structural multiplexity of the 
system. The model is illustrated in Fig. [l] 

The multiplex voter dynamics consists of a sequence of 
time steps. During a time step we perform N updates, 
and each update consists of three elementary operations, 
as follows: i) a layer a is selected at random, with uni¬ 
form probability; ii) one of the nodes i on layer a is 
chosen at random and its state s is updated according 
to the classical voter model dynamics, that is, be¬ 
comes the same as that of a randomly chosen neighbour 
of i on layer <a; Hi) if the updated node i participates in 
both layers, then the state of the corresponding node in 
the other layer /3 changes as well by setting s\^ = 

With this third operation, state changes can propagate 
across layers: indeed the presence of a fraction of nodes 
existing in both layers intertwines the voter dynamics on 
the two layers, so that in general the evolution of the 
overall multiplex dynamics might differ from the one we 
would observe on two independent networks of the same 
size. As a limiting case, the dynamics reduces to that 
of a classical voter model on a single-layer network only 
when q — 1. 

It is well known that in connected finite-size single¬ 
layer networks the voter model dynamics always reaches 
an absorbing state, where all the nodes have exactly the 
same state, in a survival time that scales with the sys¬ 
tem size N HBffiaESEI]. In networks of high effec¬ 
tive dimensionality (including random networks), when 
N oo the dynamics sustains an active disordered state 
in which nodes continue to change their state [T9] . Such 
active state is the one observed asymptotically in large 
systems, before finite-size fluctuations pull the system 


towards the absorbing state. The classical order param¬ 
eter to measure the activity in the voter model is the 
so-called interface density p{t ), defined as the fraction of 
active edges of the network, i.e. of those edges whose 
endpoints have different states. In the multiplex model 
we can define: 

?'"'«) = e> 

i j<i 

for each layer cq a = 1,2, where we denote by K ^ = 
\ Sz j the total number of edges at layer a. Notice 

that p = 0 if and only if all the nodes have the same 
state, while larger values of p are associated to active 
configurations. A second quantity of interest for a finite 
system is the average time (T) to reach an absorbing 
state of consensus. Such a time can be defined as (T) = 
J 0 °° P s (t)dt from the survival probability P s (t), i.e. the 
probability for a system to be active at time t [34] . 

For single-layer uncorrelated networks at sufficiently 
large times t > TV, P s (t) decreases exponentially P s (t) ~ 
e _2t / r [35]. Hence, (T) ~ r/2, so that the dependencies 
of (T) on system size and on the moments of the de¬ 
gree distribution are given precisely by the correspond¬ 
ing dependencies of the characteristic time r. For uncor¬ 
related networks one can find exact expressions for the 
value of the average interface density in the thermody¬ 
namic limit p smgle , as well as for the characteristic time 
to reach the absorbing state r smgle [231 ESI. An impor¬ 
tant feature is that these analytical predictions of p smgle 
and r smgle depend only on the values of the first two mo¬ 
ments of the degree distribution P(k) of the network, i.e. 
on p 1 = (k) = ^2 k kP(k ) and p 2 = {k 2 ) = ^2 k k 2 P(k ), 
and not on any other microscopic property of the net¬ 
work, and read: 


single _ /H ^ 

P ~ 3(/ii — 1) ’ 


(3) 


where p smgle is the average over surviving runs, and 


^single = (Ml ~ 1 )jAN_ 

(Ml - 2)^2 


(4) 


In the following we will focus on the values p(q,ui) and 
r(g,cj) of interface density and characteristic time of the 
multiplex voter model as a function of the average edge 
overlap of the system uj and of the fraction q of nodes 
present in both layers. 


III. NUMERICAL RESULTS 

We studied the voter model on a duplex network 
made by two random regular graphs with N nodes, each 
node having a degree equal to p. We therefore have 
Pi = (fcW) = (k^) = /q and p 2 = p 2 . In our simu¬ 
lations we fixed N = 1000 and p = 4, and we studied 
the dynamics of the system by varying the value of the 
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(a)Interface density p 


(b)Survival Probability P s (t) (c)Rescaled characteristic time r/N 


FIG. 2: (Color online) (a) Asymptotic value of the interface density averaged over the surviving runs; (b) Survival probability 
of the multiplex for different q values, for u = 0 (maroon, larger and darker markers), and u = 1 (green, smaller and lighter 
markers). The line denoted as analytics is P s m (t, q = 0, u; = 0) whose expression is given in the main text; (c) Loga rithm of the 
average characteristic time r of the approach to the absorbing state, rescaled by (1/iV) (legend same as in hg. |2(a) ). Quantities 
are functions of multiplexity g, and edge overlap u. We have considered an ensemble of 10 5 random initializations on a fixed 
duplex whose layers are random regular graphs, each with N = 1000 and (k) = 4. 


edge overlap uj and the fraction q of nodes participating 
in both layers. Since the system is finite, it will eventu¬ 
ally converge to the absorbing frozen state corresponding 
to an interface density equal to zero. Consequently, for 
each time t, we will evaluate the average value of the in¬ 
terface density only over the surviving runs, i.e. on those 
realisations of the dynamics which are still active at time 
t. We associate the asymptotic value of this numerically- 
obtained quantity with the value at the thermodynamic 
limit, given by Eq. |3|, and refer to it as interface den¬ 
sity, while being clear whether we are referring to the 
numerical or analytical results. In the case of a network 
with two layers, we have of course two values of interface 
density, one for each layer. If ( k M) = ( k the ensemble 
averages on the two layers will be equal, and hence the 
activity of an arbitrary layer is representative of the typ¬ 
ical activity of the entire system. We therefore use the 
asymptotic value of the interface density of an arbitrary 
layer to reflect the activity of the multiplex, and dispense 
with the a superscript. 

In Fig. [2( a)] we report the values of interface density 
p(q,uj) as a function of g, and for different values of uj. 
Let us consider first the case uj = 0 of no overlap be¬ 
tween the structure of the two layers of the network. 
When q = 0, i.e. when inter-layer state passing is not 
allowed, the system effectively corresponds to two identi¬ 
cal but independent single-layer voter dynamics, so that 
p(g, cJ) is in accordance with the classical analytical pre¬ 
dictions for single-layer networks (since pi = 4, Eq. © 
gives the activity p(0,0) = 0.22). On the other hand, 
when q = 1, i.e. when all the nodes participate in both 
layers and each edge exists only in one layer, the sys¬ 
tem is in all respects identical to a single-layer network 
with pi = 2 p (for p± = 8, Eq. © gives the activity 
p(l, 0) = 0.286). For intermediate values of q the dynam¬ 
ics interpolates monotonically between the two extreme 
cases, i.e. two dynamically indistinguishable voter mod¬ 
els on single-layer networks with pi m p (g m 0), and 


one voter model on a single-layer network with pi = 2p 
(q = 1). The picture changes completely as soon as co is 
large enough. In general, when the overlap is above some 
limit, then p(g, uj) is a non-monotonic concave function 
of g, with a maximum at a given value of q in [0,1] which 
depends on uj. Notice that, when uj = 1 and q = 1, i.e. if 
the two layers are identical and all the nodes participate 
in both layers, the dynamics is identical to a voter model 
on a single-layer network with pi = p. In fact, since all 
the edges exist on both layers, a node participating in 
both layers will have only p distinct neighbours, and will 
be connected to each of them on both layers. Hence, in a 
multiplex network with uj = 1, the interface density p of 
the voter model takes the same value 0.22 at q = 0 and 
at q = 1 , while for intermediate values of q the interface 
density is higher than that of a voter model on each of 
the two layers. 


We now consider the characteristic time r of the mul¬ 
tiplex dynamics, where by characteristic time we un¬ 
derstand twice the inverse exponent of the multiplex 
survival probability P™{t). We consider the multiplex 
as active if at least one of the layers is active, and 
find that for all q > 0 the survival probability of the 
multiplex does decay exponentially with some exponent 
r — T(g,u;) (%• |2(b)| ). The only exception is the q = 0 
case of the fully-disconnected multiplex, where the sur¬ 
vival probabilities of the layers are independent. In 
this case the survival probability of the multiplex is 
given by the probability that at least one of the lay¬ 
ers is active, P s m (£,g = 0 ,cj = 0) ^ P s (t)(2 — P s (t )), 
or P s m (£, 0,0) ~ 2e -2t / rSmsle — e - 4 P rSinsle . This means 
that for q = 0 the multiplex survival probability does 
not scale exp onent ially , and hence r(q = 0) is not well- 
defined. Fig. 2(c) | shows the characteristic times r(q,uj) 
obtained (the value at q = 0 is that of an approximate 
exponential fit). The trend shows a peak at small g, fol¬ 
lowed by an exponential decay with increasing number 
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of interlayer connections. The value of edge overlap uj 
controls the rate of decrease. We stress that the peak 
is not a consequence of the definition of r for the mul¬ 
tiplex, and the consequent rogue value of r at q = 0: 
the limit is truly singular, with plots for the actual av¬ 
erage time until absorption displaying the same features, 
as do the plots for the characteristic exponent layer by 
layer. We also note that the peak is robust with respect 
to system size. Therefore the slowest finite-size multi¬ 
plexes are ones where the layers are interconnected by 
the smallest number of links. Multiplexes with more in¬ 
terlayer connections will stop being active faster, as will 
wholly disconnected systems. This insight can be under¬ 
stood by realizing that in a multiplex with few interlayer 
links, the layers will most of the time function as com¬ 
pletely disconnected networks that half the time will try 
to settle into the absorbing states of consensus different 
from one another. Since this situation is now prohib¬ 
ited by the few interlayer links, the multiplex will freeze 
only when one of the layers switches over and both layers 
reach the same consensus. It is this switching behaviour 
that is responsible for the peak in r for small q. We also 
notice that this non-linear effect does not depend on the 
edge overlap uj. Here, unlike with the behaviour of in¬ 
terface density, the overlap does not change qualitatively 
the behaviour of r for increasing q. In fact, the more 
interconnected the multiplex, the higher the role of edge 
overlap: tuning up uj causes a decrease in the interface 
density and results in a longer-lived multiplex, and the 
effect becomes more pronounced with q. 


IV. IRREDUCIBILITY OF THE DYNAMICS 

It was shown in [26] that the interface density of the 
multiplex voter model dynamics as a function of q and 
in absence of edge overlap can be rewritten as the inter¬ 
face density of a voter model on a single-layer network 
having qi = (k}( 1 + q) (where k is the average degree 
of each of the two original layers), under some appropri¬ 
ately rescaled time. Here it is important to note that the 
approximation of [26, treated the interlayer connections 
as inherently probabilistic: the q parameter was a prob¬ 
ability that each node’s state gets passed on to the other 
layer. It was shown that the analytics of such a system, 
in the particular case uj = 0, are equivalent to those of 
the voter model on some properly aggregated single-layer 
network. However, as Ref. [26] only considered the ther¬ 
modynamic limit, it is not obvious whether the aggregate 
displays a corresponding rescaling of characteristic time. 
Nevertheless, it should be possible to devise an aggre¬ 
gate that results from the flattening of the multiplex into 
a single-layer network, so that the resulting network will 
have a first moment of the degree distribution equal to 
the expression fi i given above. 

The main question now is whether such a reducibility 
is possible in the most general case in which the multi¬ 
plex has edge overlap uj ^ 0. This is indeed the most 



FIG. 3: Schematic illustration of two possible aggregation 
procedures that reduce a two-layer multiplex with N nodes 
on each to a single network. The multiplex contains qN inter¬ 
layer links between associated nodes (pairwise links). These 
nodes are shown in double circles. These nodes will have twice 
as many links in the aggregates, independent of the flattening 
method. Method 1 keeps 2 N nodes in the resultant aggregate, 
whereas method 2 has N(2 — q). 

interesting case, i.e. when in a multiplex there are corre¬ 
lations between the edges at the different layers 29 . [33]. 
In the following our working hypothesis will be that, if 
the multiplex can indeed be reduced to a monoplex (in 
the variables of interest), then there exists an aggregate 
graph such that the behaviour of the voter model on the 
multiplex is completely described by Eq. © and Eq. © 
evaluated on the corresponding monoplex. Since those 
equations depend just on the number of nodes and on 
the first two moments of the degree distribution of the 
resulting single-layer network, we will now consider two 
standard aggregation procedures, and derive analytically 
the values of the first and second moments of the resulting 
degree distributions. Let us remind that, by definition, 
the first moment of the degree distribution of a graph 
is equal to /ii = (k) = ^2^ = ikP(k) 1 where P(k) is the 
degree distribution, so that P{k = k) is the probability 
that the degree of a node sampled at random from the 
graph is equal to k. The probability P(k) can be also 
written as: 

P{k) = zM (5) 

where Nk is the number of nodes in the graph having 
degree equal to k and the normalization is just the total 
number of distinct nodes. 

Given a duplex formed by two random regular graphs 
with identical degree /q with edge overlap uj and where 
a fraction q of nodes participates in both layers, we can 
distinguish two classes of nodes. The nodes in the first 
class exist in only one layer, and we indicate as k s i ng i e the 
number of their neighbours, while nodes in the second 
class exist in both layers (these are the nodes having an 
inter-layer link) and we indicate their degree as fcboth- 
The former class has degree: 

^single = Ab (b) 

whereas nodes present on both layers have degree fcboth = 
2 n when uj = 0. However, if the edge overlap is not null, 
i.e. when uj > 0, then the nodes being present in both 
layers have degree fcboth ec l ua l f° : 

fcboth = 2 fjL- qilUJ = fi {2 - quj). 


(7) 
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q q 


Method 1 Method 2 



(a)Interface density p (b)Rescaled characteristic time r/N 


(c) Error of interface density 


FIG. 4: (Color online) Comparison of the theoretical prediction for (a) the interface density and (b) the characteristic time 
based on the aggregation of the duplex. Numerical results are shown for uj — 0 (stars, black), uj — 0.5 (circles, green) and 
uj — 1 (diamonds, red), where the other parameters are the same as in Fig[2] Analytics for the aggregates are computed for 
the respective values of uj for methods 1 and 2, and normalized by N. Panel (c): Absolute value of the relative error on the 
estimation of p(q,uj) using aggregation method 1 as a function of q and uj. 


In fact, the degree of a node i present in both layers is 
equal to the sum of its degrees on the two layers (2/i) mi¬ 
nus the expected number of its edges which are present 
on both layers. This number is equal to the probabil¬ 
ity that a neighbour j of i is also present on both layers 
(which is equal to g), times the probability that the edge 
(i, j) is present in both layers (which is equal to uj) multi¬ 
plied by the number of neighbours of i (i.e., fi). Hence we 
get the correction qfiuj. In the particular case in which 
uj = 1, we get /q>oth = ^(2 — q ), while for uj = 0 we recover 

^both — 2fl. 

Let us now consider the two following distinct aggre¬ 
gation procedures. They are illustrated in Fig. [3] and 
correspond to the two most standard ways to aggregate 
a duplex into a single-layer network. The two flattening 
procedures differ in the total number of nodes and also 
in the number of nodes of degree & S i ng i e and k bo th 
they produce. This in turn changes the effective system 
size of the aggregate network, and the first two moments 
of the degree distribution. In the following sections we 
compute these quantities for the two methods, and assess 
how well the aggregates fare in describing the behaviour 
of the true multiplex. 

A. Aggregation method 1 

We can obtain a single-layer representation of a mul¬ 
tiplex by putting the layers side by side with no ef¬ 
fective node overlap. Therefore 7V S i ng i e = 27V(1 — g), 
TVboth = 2A/g, and the total number of nodes in the ag¬ 
gregate is just 7V aggr = 2 N. The degrees of the nodes 
participating in just one or both layers are respectively 
equal to k sing i e and k bo th . By using Eq. we have: 

P(k = Single) = (1 - q), (8) 

and 

(9) 


The first moment yiq ggr of the aggregated graph is then 
equal to: 

oo 

^ = E* p ( fc ) 

k=1 

= k s i ng i e P(k = k s i ng l e ) H~ k b oth.P(k = ^both) 

— /il-P(/c — ^single) T /^l (2 qUj)P{k — &both) 

= (1 - g)/ii + g(2/ii - qujfii) 

= /ii(l + g-cjg 2 ), (10) 

and the second moment 

^ = {l_q)^Pqii\{2-qujf 

= 111 [1 + 3g - Aujq 2 + uj 2 q 3 ] . (11) 

Under no overlap these reduce to 

/xf gr = Ml (l + g), 

^ ggr = /i?(l + 3g). (12) 

Note that this method produces the same effective rescal¬ 
ing of the first moment /ii as given by the analytical es¬ 
timation of the interface density of the multiplex in the 
thermodynamic limit (see Ref. [26]). 

B. Aggregation method 2 

This method reduces the number of effective nodes in 
the aggregate through treating nodes present in both lay¬ 
ers as one node. So, while for the rest of the nodes we 
still have A^ S i ng i e = 2N(1 — g), the number of ‘multiplex’ 
nodes is now half as much as in method 1, N bo th = Nq. 
The total number of distinct nodes in the aggregate is 
equal to N aggr = 2N(1 — q) + Nq = N(2 — q). The de¬ 
grees of the single and ‘multiplex’ nodes are as before 
equal to & s i ng i e and k bo th respectively. Therefore 

P(k = ^single) = "7X “T" 5 


P(k = k hoth ) = q. 


( 13 ) 
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(a)Effective degree distribution (b)Scaling of characteristic time r/N 


FIG. 5: (Color online) (a) The values of /i^ ggr resulting from the inversion of the measured interface density in the simulations 
as a function of q for different values of edge overlap (symbols) and the corresponding quadratic fits (lines), (b) Scaling of r 
with q using the data shown in Fig. [ 2 ] For each cj, a corresponds to the slope of — ln(r/N ) with q computed for q > 0. Straight 
lines are the linear regression fits for each of the two N trends, with coefficients given in the legend. 


and 

P(k = fcboth) = J^y (14) 

Consequently, we have for the first moment of the degree 
distribution /i^ sgr : 


^ = Y j kP{k) = 


2N(1 — q) x iii + Nq x p i(2 — quj) 


k=l 


N (2 - q) 


= Mi 


(2 - q 2 W ) 


(2 ~q) 

and for the second moment p 2 Sgr : 

2 (q 3 u: 2 - Aq 2 u + 2q + 2) 


M2 SSr = Mi" 


(2 -q) 


(15) 


(16) 


When the layers are uncorrelated and the overlap uo = 0, 
we have 


. aggr _ 2/ii 

Ml — 


(2 - q y 


^ = 2 ^ l+q) 


(2 ~q) 


We notice that, when q = 0 ( i.e., if we have 

two non-interacting layers), we get p^ ggr = p 1 and 
^aggr _ ^2 _ ^ while for q = 1 we have /i^ ggr = 2/ii 
and /4 ggr = 4/i? = (2 Ml ) 2 . 


C. Comparing aggregated to multiplex dynamics 

We now compare the numerical results for the multi¬ 
plex to the theoretical values given by the corresponding 


equations for the aggregates. Eqns. (|3| and Q describe 
the behaviour of the monoplex. We contrast the two ag¬ 
gregates by substituting the respective effective value for 
the first and second moments and the total number of 
nodes. In other words, we take N Af aggr , fi 1 —)► /i^ ggr , 
and /i 2 —>* M 2 SSr , where the effective values differ depend¬ 
ing on the aggregation method. 


The results are reported in Fig. [4] Consider first the 
behaviour in the thermodynamic limit, described by the 
interface density p. Both methods result in aggregates 
whose qualitative behaviour with uj and q c orresp onds 
to the trend observed in the multiplex (Fig. |4(a)| ): for 
small q the system becomes more active with increasing 
number of interlayer edges, whereas after a certain point 
activity may decrease, depending on the whether or not 
the edge overlap is significant enough. However, neither 
method gives correct quantitive understanding for a gen¬ 
eral uj and q > 0. This can be seen in Fig. 4(c), which 
compares the performance of the two methods through 
computing the absolute value of the relative error on p as 
heat map of q and uj. Method 1 in general does a better 
job than method 2 at being systematically consistent, al¬ 
beit suggesting higher values (the semicircular drop in the 
error to zero, observed in method 2, comes about from 
the crossing of the analytical and the numerical trends). 
As we have seen before, the analytical results for Method 
1 (at uj = 0) correspond to a multiplex with probabilis¬ 
tic interlayer connections of intensity q. Therefore, at 
uj = 0, comparing the numerical results for the multi¬ 
plex, and the analytical trend for method 1 can be used 
to gain insight into the difference induced by an alterna¬ 
tive method of inteconnecting the layers of the multiplex 
(in fact, when viewed in this way, the small magnitude of 
the differences becomes more surprising than their pres- 


























ence). The main quantitative differences arise for a wide 
range of intermediate overlap values, for medium-to-large 
interlayer connectivity - precisely the region of parameter 
space motivated by real-world systems. The conclusion 
we draw from this is that, as long as edge overlap is 
taken into account, standard aggregates can only inform 
on some qualitative features of large multiplex systems, 
but cannot capture quantitatively the behavior of real 
multiplex networks. 

Consider now the characteristic time of th e app roach 
to absorbing states of finite-size systems (Fig. 4(b) | ) . The 
difference at q = 0 is due to r(q = 0 , cj ) not being de¬ 
fined for the multiplex; yet the aggregate timescale is 
close. The monotonic decrease for q > 0 shown by the 
numerics is captured only by method 2, with coincidence 
of r values for q == 1. In fact, method 2 does a better 
job at both qualitative and quantitative results, unlike 
method 1, which appears to work best in the thermody¬ 
namic limit (this aggregation method results in a much 
slower system for large q than the observed duplex). The 
main discrepancy, however, is that neither of the meth¬ 
ods capture the jump of r at small q. Therefore, although 
the two standard aggregation procedures can come close 
to describing the qualitative features of long-term duplex 
activity, they are not sophisticated enough to capture the 
long timescales associated with sparsely interconnected 
systems. In the following Section we will see whether 
and how it is possible to devise more complex aggrega¬ 
tion procedures in order to reproduce quantitatively the 
dynamics of the multiplex voter model. 


V. NON-LINEAR MULTIPLEX EFFECTS 

In principle, the results of the previous Section can¬ 
not absolutely rule out the possibility that there exist 
other aggregation methods reproducing the behaviour of 
p(q,uj) observed in the multiplex voter model. Our hy¬ 
pothesis is that the deviations from the theory found for 
uj > 0 are due to the additional non-linearity induced by 
inter-layer state copying made possible by the presence of 
a fraction of inter-layer edges. In order to better investi¬ 
gate these non-linear effects, we inverted Eq. (|3| to com¬ 
pute the effective value of the average degree p J1 ggr (q,uj) 
of an ideal aggregated single-layer network able to repro¬ 
duce the observed value of p(q, uj) for each value of uj and 
q. In formula: 


^ sgr (q,cv) 


3 p(g,u>) - 2 
3 p(q,u) - 1 


(17) 


As made evident by Fig.[5j for each value of edge overlap, 
/i)' gsr ((/, oj) can be fitted very well by a quadratic polyno- 
mial in q. However, the actual values of the coefficients of 
the fit depend on uj in a non-trivial way. We can formally 
write: 


p\ ggr (q,uj) = F(uj)q 2 + G(uj)q + H(uj). (18) 


The problem now is to find an expression for F(uj), G(uj) 
and H{uj). We started by making an ansatz for the func¬ 
tional dependence of those coefficients on uj, and then 
fitting these functions, for each value of q, by using sev¬ 
eral realisations of /i^ sgr (g, uj) corresponding to different 
values of uj. 

We found that for all the values of uj the quadratic fit 
of yiq ggr (g,u;) yields H(uj) = p, so we just focused on the 
other two coefficients. The values of F(uj) and G(uj) are 
reported in the panels of Fig.[6]as black circles. It is then 
reasonable to assume that F(uj), G(uj) are polynomial 
functions of uj. We found that in order to accurately 
reproduce the behaviour of /i^ ggr (g, uj) in the whole range 
of uj, both F(uj) and G(uj) should be at least third-order 
polynomials in uj, as made evident by the plots reported 
in Fig. [6j 

Notice that the predictions of /i^ ggr provided by the 
two theoretical arguments reported above, based on the 
linear superpositions of the two layers, contained only lin¬ 
ear terms in uj. However, the fit of p^ ggr (q,uj) confirms 
that, when uj > 0, the behaviour of the multiplex voter 
model is the result of a highly non-linear combination 
of the two layers, suggesting that a trivial single-layer 
equivalent of the multiplex voter model dynamics does 
not exist, especially when the underlying multiplex net¬ 
work is characterised by a non-negligible overlap. 

That simple aggregation procedures do not produce 
the observed scaling is also evident by examining the 
timescales on which finite systems approach absorbing 
states. For q > 0, the characteristic time r(q,uj) is 
an exponentially decrea sing function o f q, t (q,uj)/N ~ 
e -qa(uj,N)^ a > o (fig. 2(b)). Figure 5(b)| shows that 


a(uj,N) = — auj + b, a > 0, and that the system-size 
dependency does not enter into it. Hence, r(q,uj)/N ~ 
e q(ua-b) ' Thus increasing the edge overlap uj results in 
longer-lived systems, while adding more interlayer links 
produces the opposite effect. This was additionally con¬ 
firmed by examining the behaviour of the rescaled time 
until absorption (T), which showed the same qualitative 
trend and almost identical a and b coefficients. 


VI. DISCUSSION 

Multilayer networks allow to extend the applicability of 
network theory to more realistic contexts in which nodes 
are connected through concurrent interaction patterns 
of different kinds. However, a fundamental open ques¬ 
tion to answer is whether the added complexity yielded 
by multilayer networks is really needed to model net¬ 
work phenomena, or if instead there exist simple ways of 
representing multiplex dynamics through appropriately 
constructed processes occurring on appropriately con¬ 
structed single-layer networks. We have investigated here 
the problem of reducing the multiplex voter model to an 
equivalent single-layer dynamics. We have considered the 
predictions about the level of activity of the multiplex 
voter model in the thermodynamic limit, as measured 
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FIG. 6: The black dots in each panel represent the coefficients of the quadratic fit F(uj)q 2 + G{ui)q + 1 of fi\ sgr (g, cj) for different 
values of u. Notice that both F(uj) and G(uS) are non-linear functions of uj. The different lines correspond to fits of F(uj) and 
G(uS) using d -th order polynomials. Interestingly, a good fit is obtained only for d > 3, meaning that the presence of inter-layer 
state copying is introducing highly non-linear effects. 


by the interface density, as well as the time to reach the 
absorbing state for finite systems. We have found that 
results for the interface density based on single-layer ag¬ 
gregated graphs are accurate only when there is little or 
no interaction between the layers (q sufficiently small), 
or when there is full connectivity (q ~ 1) and the edge 
overlap uj is either 0 or 1. For the complementary broad 
range of parameters, any standard aggregation procedure 
can only give some qualitative information about the in¬ 
terface density, but fails to reproduce the quantitative 
details. We showed that edge overlap and multiplexity 
have two opposite effects on the long-term dynamics of 
the multiplex voter model, and in particular that an in¬ 
crease in the value of edge overlap can counter the action 
of increasing the fraction of interlayer links, leading to 
an overall decrease in the interface density of the mul¬ 
tiplex. In fact, we showed numerically that any equiv¬ 
alent single-layer representation of the multiplex voter 
model dynamics entails the construction of an aggregate 
network which is a highly non-linear combination of the 
original layers. 

These results will be found at the same time surprising 
and interesting by all the researchers aiming at modelling 
social interaction in real-world scenarios. As a matter 
of fact, it has been recently shown that multilayer so¬ 
cial networks are normally truly multiplex , meaning that 
they are characterised by non-negligible values of edge 
overlap [30] and by an intermediate level of multiplex¬ 
ity [29l 3T. 32^. And as we have shown in this work, the 
single-layer voter model approximation is qualitatively 


(when considering the characteristic time) and quantita¬ 
tively (when considering both the characteristic time and 
the interface density) inaccurate when the edge overlap 
uj and the degree of multiplexing q are far from their 
extreme values 0 or 1. This is reflected in the average 
timescales of consensus for finite systems, as we found 
that a multiplex with very few interlayer connections 
takes the longest time to reach consensus, much more 
than would probabilistically be needed by two discon¬ 
nected monoplexes. This nonlinear effect is not captured 
by either of the simple aggregation procedures proposed 
in the paper. 

The case of dynamical irreducibility of a multiplex pro¬ 
cess presented in this work raises the important question 
of whether other unknown phenomena might be lurking 
in the multilayer structure of real-world systems. This 
question, together with the insights about the intrin¬ 
sically multidimensional nature of the multiplex voter 
model, represents a stimulus to perform further research 
along these lines. 
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